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MICROCOPY RESOLUTION TEST CHART 
NATIONAL BUREAU OF STANDARDS 
STANDARD REFERENCE MATERIAL 1010a 
(ANSI and ISO TEST CHART No. 2) 


3. DISPLACED CAPACITY 
3484 Introducticn 


The installed capacity of a conventional power generating system should 
be such that the system is capable of meeting a certain maximum demand 
(possibly with uncertainty in the forecast) at some level of re- 
liability. This level is expressed in a reliability index which can be 
determined from statistical considerations of load demand and available 
power. With the integration of wind power, with its stochastic nature, 
the question may be posed how much conventional capacity may be removed 
from the system in order to maintain the required level of re- 
liability. This displaced capacity is also called capacity credit. In- 
stead of removing conventional capacity, one can detc<rmine the in- 
crement of maximum demand which could be served with maintained re- 
liability. The capacity credit of wind power is then equal to the 
amount of conventional capacity which could serve the same increment of 
maximum demand. These two interpretations of capacity credit - both 
considering the conventional equivalent of wind power as regards 


reliability - yield almost the same results 7]. 


Three reliability criteria will be considered in this report (which, 
however, are not independent of each other for a given demand and 
supply pattern): 

- The probability that the available power during annual peak load is 
insufficient to serve this load should not exceed a predetermined 
value (say 1%). This probability is called loss-of-capacity probability. 
The same criterion can be applied for periods shorter than a year, 
such that scheduled maintenance of power plants can be taken into 
account adequately. 

The expected total time fraction with insufficient available power 
should not exceed a predetermined value. This expected fraction is 
called loss-of-load probability and can be expressed in hours per 
year. 

The displaced conventional capacity, calculated with the above criteria, 

is the subject of sections 3.2, 3.3 and 3.4 respectively. The effects 

of various assumptions will be demonstrated for wind plants at Den 

Helder only. The effects of dispersed siting are dealt with in section 


). 2. 





3.2. The loss-of-capacity method for the annual peal 


Consider the probability that a power system consist 
number of plants has a certain capacity available at 
Each plant of capacity C; is either available, with 
not available (due to break-down), with probability c 
general q, is օք the order of 102). The probabilitie 
are mutually independent. For a system consisting of 
large number of plants (at least more than 20 units) 
density function of the available capacity p(C) may 
a normal or Gaussian distribution and therefore cha 
parameters only: 

- mean available capacity: 


- varianc ¿č available capacity: 


We may now require that the probability of the avai 
being less than peak demand — does not exceed a 
12). Thus: 


max 
f  p(O dC = a. 
0 


In the normal distribution approximation this means 


of the system should be chosen such that: 


— = ローC(o) . 0 。 


where c(a) is the value above which the area under t 
unit variance is equal to a (a-quantile); e.g., c(0. 
c(0.03) = 1.88. 

For the Dutch system that will exist in 1985, as giv 


unit forced outage rates given in 115 | , the followin 


= 14037 MW; o = 690 MW. 


Having found these values for the conventional syste 
calculate the permissible increment of Lax when win 


capacity is incorporated. The probability distributi 


. load 


ing of a great 

a certain moment. 
probability Pi, Or 
լ Հ 1-pi (in 

Տ of the plants 

|. a sufficiently 

չ the probability 

be approximated by 


acterised by two 


[3.8.1] 
(3.2.2) 


able capacity 


mall value a (say 


(3.2.3) 


that the capacity 


(3.2.4) 


he Gaussian with 


01) = 2.33, and 


en in |17|, with 


g 1s found: 


m, one can 
d power of specific 


on function of 


the combined system, p(C), is now obtained by convolution of the 
normal distribution function of the available conventional power with 
the probability distribution function of the wind power. The latter 
distribution function is the wind power duration curve as described 
in section 2.2. With eq. (3.2.3) the new Lmax For the same value of 

a can be found. The increment օք Lmax ԼՏ called the load carrying 
capability (LCC) of the added plant. 

These calculations can be repeated for the addition of one or more 
conventional plants - with a binary probability distribution for a 
single plant and a binomial distribution if more plants are involved. 
The wind power capacity credit is equal to the conventional capacity 
which has the same LCC as the wind plant. 

Fig. 7 shows LCC values for wind plants and conventional units. The 
LCC values are determined by the properties of the added plant, the 
standard deviation o of the existing conventional system and the 
reliability index a. For very small values of added capacity (wind 

or conventional) the LCC equals average available power. For very 
large values of added capacity LCC tends to a limiting value; for this 
limit the time fraction with zero available power is the only 


property of interest of the added plant (wind or conventional). 


In order to obtain more accurate figures for the required conventional 
system reserve, KEMA |15| has introduced some refinements in the method 
discussed in the previous section. The same improvements can be 

adopted for the investigation of wind power integration. In this section 
we will consider the effects of the uncertainty in forecast maximum 
load, of the interconnection of the domestic grid with those of 


neighbouring countries, and of time-dependent unit availabilities. 


The effects of scheduled maintenance will be discussed in section 3.3. 


In their expansion planning studies |15| the Dutch utilities take 
into account an uncertainty in the forecast maximum load կաւ" This 
uncertainty is assumed to have a normal distribution with standard 


deviation equal to 3.27 of L ax? i.e. 367 MW for the year 1985. 


Thus both demand and supply have normal distributions and this is also 
true for their difference. The standard deviation c to be used in, 
e.g., eq. (3.2.4) then increases for the year 1985 from 690 MW to 

782 MW. 

The effect of the uncertainty in the peak load on the calculated 
capacity credit of wind power is shown in fig. 8. Here it was assumed 
that the load carrying capability of conventional units to be replaced 
by wind power was simply 902 of their capacity, which is a fairly 
optimistic value (see section 3.2.1 and fig. 7); the resulting capacity 
credit will, therefore, be rather pessimistic. In fig. 8 curve A 
represents the capacity credit as a function of the installed wind 
power (in MW) with a = 0.01 and o = 690 MW and thus without forecast 
uncertainty. Curve B replaces curve A if this uncertainty is taken 


into account. 





Grid interconnections 


For the evaluation of the interconnection of the domestic grid with 
those of neighbouring countries (Germany & Belgium) the same philosophy 
as in |15| is adopted. In this philosophy, the foreign supply system 
is assumed to be identical to the domestic system. Both systems have 
to serve identical loads. All power plant outages are supposed to be 
uncorrelated; the two loads are assumed to be fully correlated. 
Capacities of tie lines connecting the two grids do not pose any 
power transfer limitations. If one requires a reliability level 

a = 0.01 for the joint system, the effect of the interconnection can 
be interpreted as a relaxation of the required reliability level for 
the domestic grid, which can be determined as follows. 


For the domestic system eq. (3.2.4) reads 





L = ローc(o) .o 。 c = V690? + 367? , (3.2.5) 
max 
where 690 (MW) is the standard deviation of available power and 367 (MW) 
the forecast uncertainty. 


For the joint system eq. (3.2.4) reads 





2L - 2u - c(0.01) . c., g. = V5 . 6902 + 27,3672 . (3.2.6) 
max J J 


Combination of the two equations yields 


c(a) = 1.82 + a = 0.035. 


In fig. 8 curve C replaces curve B if the domestic loss-of-capacity 
probability is allowed to be 3.52. 

It should be noted that curve C is based on the assumption that, if 
wind plants would be installed in the other countries, there would be 
no time correlation between the power outputs of the wind plants in 


different countries. 





In reliability studies the Dutch utilities make use of unit avail- 
abilities appropriate to hours of high demand 15| . In gathering 
statistical information on these availabilities they do not consider 
availabilities during night hours. The same can be done with respect 

to wind power availability, in which case the seasonal and diurnal 
variations are of importance (see fig. 3). Peak load in 1985 is supposed 
to occur in December. We therefore determined the wind power duration 
curve for the hours between 6.40 a.m. and 10.40 p.m. ir December. 


With these wind data curve C in fig. 8 is replaced by curve D. 


All refinements discussed up to now appear to have a smali positive 
effect on the capacity credit of wind power. Nevertheless their sum 


is significant. 


There is some danger in using wind power distribution functions based 
on long time series of wind data. In fig. 3 one can see a fluctuation 
in average annual available wind power. The effect of this variation 
on capacity credit is shown in fig. 9, where the "average" CCR is 
compared with the CCR's based on wind data of the years 1963 and 1966, 
being the poorest and richest aeolian years, respectively, of the whole 
period of observation. One may wonder what the loss-of-capacity 
probability would be in such extreme wind years if a capacity dis- 
placement would have been realized based on "average" wind availability. 
Table 3 shows these values for two wind plants with installed 
capacities of 1000 MW and 3000 MW, respectively. Even for 3000 MW wind 


power the effect is quite small, but the question whether these annual 


Ա 


variations in reliability level are acceptable goes beyond the scope 


of this study. 


Table 3. Capacity credit CCR of wind plants based on various wind data 





of Den Helder, and loss-of-load probability a for wind data 
of extreme years. Load forecast uncertainty and grid inter- 


connection have been taken into account. 





a (ん) for CCR based 


CCR (MW) for a= 17 on all wind data 





installed wind power 1000 MW | 3000 MW 1000 MW | 3000 MW 





all wind data 258 
1963 wind data 215 
1966 wind data 292 























3.3. The loss-of-capacity method for 4-weekly peak loads 





The seasonal load variation is rather small in the Netherlands. There- 
fore, scheduled maintenance cannot be confined to periods with low 
demand only, without adversely affecting system reliability during 
these periods. Dutch utilities plan their maintenance scheme in such 

a way that system reliability remains about the same over the year. 

In this planning procedure the year is divided into 13 periods of 4 
weeks. For each 4-weekly peak load the loss-of-capacity probability 
should be about 1Z. On the average, thermal units are on scheduled 
maintenance during 8.52 of the time. Owing to this variation of con- 
ventional capacity available for productior the load carrying capa- 
bility LCC of wind power will also vary somewhat. Furthermore, if LCC 
values per period are to be determined, it is reasonable to take into 
account the seasonal variation in available wind power. We therefore 
determined wind power distribution functions per period for day-hours 
(6.50 a.m. - 10.40 p.m.). Using these functions, together with the 
maintenance scheme for 1985 as given in |15|, LCC values per period can 
be derived. Table 4 gives an impression of the variation of LCC over 


the year. The 13 peak loads and their mean value (Ly), as well as the 
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average of the 13 LCC values (LCC) are also shown. 
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Table 4. Load carrying capabilities LCC (MW) of two wind plants 





(installed capacity C sat? for 4-weekly periods. Loss-of- 


capacity probability a - 


0.01 


(including the effects of 


forecast uncertainty and grid interconnection). Wind data of 


Den Helder (period-dependent, day hours only). 





























period | ww | inst - 
M^ 1000 3000 
| 11412 24] 388 
2 | 11374 | 276 429 
3 | 1145] | 298 476 
4 | 11103 ] 279 433 | 
5 | 10982 | 226 364 | 
6 | 10629 | 210 361 | 
7 | 10345 | 248 402 | 
8 | — 10141 | 205 344 | 
E | 9651 | 244 387 | 
10 10465 | 188 300 | 
11 10629 229 351 | 
12 11172 | 257 394 | 
13 | 11472 | 245 | 377 | 
| | | թ 
— | 10833 | 242 | 385 
! value | | 


In order to obtain a value for the capacity credit of wind power we 


make the following simplifying assumptions: 


- The maintenance scheme for conventional units will not be adapted 


to the seasonal wind power variation. In that case it can be shown 


that, if - with wind integration - all expected load values in the 


year under consideration would increase by a fraction equal to 


LCC/Ly , the average reliability would remain on almost the same 


level (although the reliability index per period would be altered 


gomewhat, see below). 


- As before, it is assumed that conventional capacity to be displaced 
by wind power has a load carrying capability equal to 902 of in- 
stalled capacity. It is further assumed that the maintenance scheme 
for this capacity could have been optimized in such a way that in 


each period the relative load increase to be served by this capacity 


A. 


power generating systems. 


is equal to the above mentioned value LCC/Ly. 


With these assumptions, and using an average maintenance time fraction 


of 0.085, it is easily shown that the conventional capacity that can 


Table 5. Loss-of-capacity probabilities a. of three 





be displaced is equal to 


CCR = LCC 





1000 MW 
wind power 
a (2) 


3000 MW 


wind power | 


a (Z) 





0.9 x(1-0.085) ' (3.3.1) 


The assumptions made here with respect to the conventional units will 


again lead to pessimistic values for the wind power capacity credit. 
If the effects on CCR of scheduled maintenance are taken into account 


in this way, curve D in fig. 8 is replaced by curve E. 


It is still necessary to show that the reliability levels Գ. for each 
period remain acceptable. In table 5 a comparison is made between a; 
values for the 1985 supply system without wind 115] and systems with 


1000 MW ard 3000 MW wind power if capacity displacements according to 


mo Ww OO Վ Հ un と を vv 9 


eq. (3.3.1) are realized. One can see that the fluctuation of a is a 
bit more pronounced after wind power integration. The reason for this 


effect being so small is that even for large penetration levels of wind 





power the capacity displacement is only a very small fraction of total | 


mm 


l. 
1.16 
(.0004) 
1.40 
0.95 
1.31 
0.65 





0.13 
0.89 
0.84 
1.12 





0.12 
0.88 
0.77 











installed capacity. A small adjustment of the maintenance scheme could | 
geometric 


further improve overall system reliability. All a values given in table 5 


mean x | 
L 


meet with the three reliability criteria mentioned in | 15|: 








l. The geometric mean of a should be less or equal 12 

d. o. should not exceed 27 in any period. 

3. գ. should not exceed 17 in periods with relatively high demand and, 
therefore, relatively little scheduled maintenance (these are 


periods 1,2,3 and 13). 





In taking this average period 9 must be omitted. This period co- 
incides with summer holidays. A low power demand and almost no 
plants withdrawn for maintenance cause the loss-of-capacity pro- 
bability to drop to an extremely low value. 


Finally, it should be mentioned that maintenance and technical outages 
əf wind turbines have not been considered yet. The main reason is 

that statistical information on these subjects is lacking. However, 

if wind turbines have an average outage time fraction equal to զ 
(planned plus unforeseen outages), it can be argued that fig. 8 still 
contains correct normalized capacity credit values if Cinst in this 


figure is interpreted as the "effective" installec capacity 


Ceff = (1-9) .Cinst: 


3.4. The loss-of-load method 





In contrast to the loss-of-capacity method the loss-of-load method 
considers the probability to be in a loss-of-load situation at any 
moment of the period under investigation. It is a generalization of 
the well-known procedure in which, for a known load duration curve 
(LDC), the plants are set on-load one by one (or in groups) in order 
to fill the area under the LDC, starting with base load and ending with 
peak load. In the generalized method acccunt is taken of the fact 
that units will have availabilities of less than 1002, by using the 
LDC of the "effective load". The effective load is defined as the 

sum of actual load and capacity on forced outage. In the same way as 
described above the units are set on-1oad one by one, and with the 
last unit in operation the remaining area below the effective LDC is 
a measure of the demand not met; the loss-of-load probability (LOLP) 
as well as the curtailed energy (loss-of-energy) is then eesily 
determined. In comparison with a deterministic method this procedure 
also has the advantage that results are obtained concerning capacity 
and load of peaking units (stalled reserve). 

The theory can be found in literature |16| and is often given for two- 
state systems (a unit is either available or on outage). However, the 
generalization to include any arbitrary probability distribution is 
easily realized and thus allows investigation of the integration of 
wind power,into the utility grid. 


An abstract of the theory is given in the Appendix. 


The loss-of-load method can be used for cost optimization calculations, 
by appropriate ordering of the plants. For the calculation of the re- 


liability index LOLP the loading order of the plants is irrelevant. 





Chapter 4 will dea 


attention to the o 


Results of wind po 
the loss-of-load p 
sults based on the 
The two sets of re 
installed wind pow 
terion 1s a bit hi 
generating power a 
peak hours。 which 

of this on the loa 
wind plants than f 


power. 


3.5. Effect of dis 





Up to now wind dat 
parison is made of 
locations mentione 
Texel and Goeree t 
steeply with incre 
for Den Helder and 
availability above 
this effect will b 
Here use was made 

cribed in section 

normalized to the 

Even for the four ` 
relatively close t 


on the capacity cr 


3.6. Conclusions 





The concept of dis 
an existing conven 
meaningful provide 
changes in the com 


in this chapter yi 
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1 with the cost optimization calculation and special 


rdering of the plants will be given there. 


wer capacity credit calculations based on maintaining 
robability are given in fig. 10, together with re- 

loss-of-capacity method for the annual peak load. 
sults compare very weil. For intermediate values cf 
er the capacity credit according to the LOLP cri- 
gher, which may be explained as follows. Available 
lso adds somewhat to system r2liability during off- 
is taken into account in the LOLP method. The effect 
d carrying capability is more important for large 


or small conventional units to be displaced by wind 


persed siting 





a of Den Helder only were used. In fig. 11 a com- 
the capacity credits of wind plants at the four 
d in section 2.2. For the North Sea locations 
he relative capacity credit CCR/Cinst drops less 
asing wind power penetration than these functions 
Urk do, which is due to the more regular wind 
sea. If dispersed siting of wind plants is applied 
e still more pronounced, as is shown in fig. 12. 
of the combined wind power duration functions des- 
2.2. For ease of comparison the CCR values were 
average power output of the (combined) plants 
weather stations considered here, which are located 
o each other, dispersed siting has a positive effect 


edit which should not be neglected. 


placed capacity, when wind power is integrated into 
tional electricity generating system, seems to be 
d the integration does not bring forth significant 
ventional system. The different methods discussed 


eld fairlv consistent results with respect to the 


displaced capacity. 

The variations of wind power availability from year to year, for the 
period investigated, are mild and do not substantially affect the 
capacity credit to be assigned to wind plants. There is a small but 
positive correlation between load demand and wind availability. Dis- 
persed siting of wind plants has a clear favourable influence on their 


capacity credit. 


Further study is required with respect to the following problem areas: 

- The extrapolation to turbine hub height of wind speeds measured 
at ground level should be made more reliable. 

- The interconnection of the domestic grid with foreign ones was 
taken into account in a simple and rather arbitrary way. A further 
evaluation of the effect of this interconnection seems to be ne- 
cessary. 

- Energy storage has not been considered. The LOLP method (section 


3.4) seems to be suitable to investigate the effects of storage 


facilities on system reliability |18|. One conclusion that can be 


drawn from literature is that all system components (storage * wind 
power * conventional plant) should be considered together since their 
contributions to system reliability are mutually dependent. 

System reliability has been expressed in a single numerical value. 
Several topics such as load following capability of conventional 
units and wind power predictability have not been considered here, but 
may affect system reliability. Meanwhile, methods which give a more 
complete description of system reliability have been applied by us 


in studies to be reported elsewhere. 
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They were all first scheduled to operate at their technical minimum Comparing the tvo strategies, one can observe a signi 


Table 6. Simplified scheme of the power system in 1985. 





4. WHA IS SAVED BY WIND ENERCY AND HOW TO ASSESS IT (30-40% of rated power) to serve low system demand, and next allowed for the oil/gas fired plants. In variant 2 the produc 





smaller and less efficient plants is much larger than 





to increase their output for load-following. Finally, under normal 

4.1. The production of individual plants i Size | Number | Installed Forced Maintenance 

Fuel in of capacity | outage rate | unavailability 

In the previous chapter tLe effect of wind power integration on total | MW | Plants in MW in ん | in weeks/year 
| | 


uranium | 12 


There is also a shift in fuel savings due to wind pow 





load-following conditions the fossil fired units were not allowed to 





operate at their rated capacity, but only at about 907 of this value; In variant 2 the savings are partially shifted from t 





required conventional capacity was considered. Another consequence of the last 107 was used only for system peak demand and in emergency to the larger ones. Furthermore, in this strategy abo 








situations. In this way a reasonable modelling of spinning reserve able wind energy of the 4000 MW wind power plant cann 


ov Oc 


tiis integration is, that the load factors of conventional plants | | 
uranium | 12 


coal | 10 





wili change. This chapter will deal with the production of the indivi- could be established. The loading order and programmed operation of and is wasted. 


dual plants. 0 は は 7 | 10 the plants are shown in table 7. Results for this strategy are given Although some operation restrictions have been taken 
The calculation method applied in this chapter is the LOLP method men- oil/gas | 12 in fig. 13, variant 2. variant 2, which lei-to less attractive figures for t 


tioned in section 3.4. and explained in the Appendix. As is indicated to wind power, the results may still be too optimisti 


oil/gas 12 


Table 7. Load ordering of the conventional power plants and wind to an increased number of stops and starts and power 


in the Appendix, the power duration curve of the individual plants and 
potential (variant 2). 





oil/gas 12 


the energy produced can be determined with this method, as soon as balance the fluctuating and partly unpredictable wind 


oil/gas 12 





the loading order of the plants is chosen. been dealt with. The time-dependent behaviour of the 


մ pd 


Size umber of | Load in percentage 
ւո MW plants of rated power 


oil/gas 12 
It is therefore possible to investigate the change of the total fuel oil/gaà 12 can be modelled in more detail with the so-cclled fre 


consumption and the fuel mix with varying strategies of the loading order. 


In all calculations to be discussed in this chapter use was made օք 
the forecast electricity demand in 1985 shown in fig. 6. Wind plants 
were a:sumed to be located in Den Helder and Urk; the wind power dura- 
tion curve for the combined output of these two plants is given in 


Fil. 2. 


In this section some results are given for a conventional generating 
system which approaches the planned system for 1985. This simplified 
system is shown iu table 6. The simplification mainly concerns the 
number of capacity classes of the bi-valent gas/oil fired plants. 
Plants are distinguished by age, under 15 years and between 15 and 20 
years oí age; plants older than 20 years were left out. 

In the calculations the unit forced outage rates were augmented to 
include the outages for scheduled maintenance. These outage rates were 


| 


taken from ref. |15|. In order to investigate the influence of the 
chosen unit scheduling strategy, two loading orders were considered. 
First, a very simple scheme was assumed. All available wind power was 
accepted first. Then the conventional units were committec in the 
order shown in table 6, 1.6. the order of increasing pro- 

duction costs. All units were allowed to operate at once at their 
rated capacity when necessary. Calculated results for this strategy, 


and for varying penetrations of wind power, are shown in fig. 13, 


variant |. 
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oil/gas 25 12 























- 
» 


a) Հ 15 years 


b) » 15 years, Հ 20 years. 


The figure shows that the total capacity of the nuclear and coal fired 
plants is so small that their production is not curtailed by even 
4000 MW installed wind capacity. Fuel savings are obtained for oil/gas 
fired plants and in particular for the smaller and older units with 
worst fuel economy. One has to bear in mind, however, that the calcu- 
lated fuel savings will be too optimistic since no account has been 
taken of any possible operating constraint or load-following penalty. 
In variant 2 a more elaborate strategy of plant operation was assumed. 
First, nuclear plart was considered inflexible, operating at rated 
power whenever available. Secondly, since start and stop procedures 
for the large fossil fired plants are undesirable from a technical and 
economic point of view, these units (» 375 MW) were kept on-line even 
if available wind power would have allowed the shut-down of some of 
these units. 

Further,a programmed operation of the gas/oil fired plants was applied, 


in which the power levels of these units were changed stepwise: 





uranium 
uranium 
coal 

o1l/gas 
wind 

Oo11/gas 
o11/gas 
oil/gas 
o11/gas 
o11/gas 
coal 

o1!/gas 


oil/gas 





i 


450 

50 

600 
575-375 





| variable 
| 200-125 a) 
200-125 Ե) 
575-125 Յ) 
200-125 Ե) 
25 a) 
600 
575-125 a) 


| 200-25 Ե) 
| 











100 
100 
90 








a) plants not older than 15 years. 


b) plants older than 15 years. 





method, a somewhat more elaborate probabilistic metho 
progress to use this method in wind integration studi 


these investigations will be published elsewhere. 


Further application of the LOLP method may also inclu 
co-generation plants for district heating and industr 
these plants the generated power is not necessarily d 
electricity demand (in the above calculations the co- 
for district heating have not been taken into account 
bution in 1985 will be negligibly small), Also the mer 


storage can be investigated with this method (see sec 


In the method discussed the electricity demand and wi 
Supposed to be uncorrelated. This may proof unsatisfa 
periods, one year for example, are considered. Howeve 
any correlations can be quantified by performing calc 
the probability distribution of the momentary differe 
demand and wind potential |5|. An alternative way is 

year in shorter periods, such that the fluctuations w 
the mean load of a period and the fluctuations of the 
with respect to the mean of the same period are stati 
dent. Shorter periods are anyhow necessary if accurat 
required, e.g. if maintenance schemes are to be taken 


detail. 
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4.2. System optimization 





In the previous section the LOLP method was used for a given conven- 
tional system with predetermined load ordering scheme. A realistic 
ordering scheme should yield the lowest production costs that are 
possible within the operating constraints. 

For long term planning it will be useful also to optimize the plant 
mix of the supply system, with or without wind power. In that case it 
might be interesting to optimize the plant mix in such a way that 
total costs are minimized (even if such a criterion need not be very 
realistic, as will be shown later). In this section results will be 
given of a few rather simple optimization calculations, which may show 
t - general trend that can be anticipated if wind power is to be in- 
tegrated. 

For these calculations a limited number of plant types was specified, 
distinguished by fuel, rated power, non-availability, and fixed and 
variable costs as shown in table 8. 

The cost figures in this table are only -easonable guesses; their ab- 
solute values are not important for the present calculations. In the 
process of composing the optimal plant mix, operating constraints were 
neglected. Thus, for a given demand, wind power was accepted first, 
followed by the conventional plants as explained in section 3.4. and 
in the Appendix. After each step the plant type of the unit to be se- 


lected next should have lowest marginal costs: 
CF/t * CV. (4.2.1) 


CF is the annual fixed costs of the plant per MW. 
CV is the variable costs per produced MWy. 


t is the time fraction with remaining demand > 0. 





Table 8. Properties of five conventional power plant types. 





Size Percentage beoe: Variable 


i costs 
ւո of non- 


costs 
MW availability Ct/ kW ct/kWh 
per hr 





uranium 10 6. .0 
coal 10 4. 
oil/gas 10 3. 
oil/gas 10 
oil/gas 























In the probabilistic method used here two refinements were incorporated 


in the expression for the marginal costs: 

- The fixed costs part (=) was increased to account for unit non- 
availability and was replaced by ՀՈ»: զ being the non-availability 
probability. 

It should be taken into account that the different unit types do 
not equally contribute to the system load carrying capobility and, 
therefore, do not equally reduce the required peak load capacity. 
This effect was accounted for by reducing the marginal costs with 
an amount equal to the fixed costs of the displaced peaking plant 
capacity (peaking plants are the units with lowest fixed costs), 
again properly normalized to 1 MWy. For this correction the load 
carrying capabilities of the different plants, calculated according 


to the loss-of-capacity method (see section 3.2. and fig. 7), were 


approximated by an analytical expression 15,19]. 


The process of adding new units was terminated as soon as the remaining 


loss-of-load probability was less than a predetermined value. 


For the forecast electricity demand of 1985, as shown in fig. 6, the 

following situations were considered in the calculations: 

a) An optimized generating system without wind. 

b) Integration of 4000 MW wind power with the system found under a) 
(in the same way as variant ]l of section 4.1.). 

c) An optimized generating system consisting of 4000 MW wind and a 


mix of conventional plants. 


The main results concerning capacities, productions and costs of the 
conventional plants are given in table 9. Costs and production figures 
in this table are given per hour and should be interpreted as annual 
average values. Comparison of the results for case a) and case b) 
shows that integration of wind power into an existing conventional 
system may effect savings for the higher priced fuels (if no operating 
constraints would be active): the average variable costs per kWh have 
dropped somewhat. This is in agreement with the results of section 
4.1. However, total costs per kWh are higher in case b) because the 
fixed costs per kWh have increased substantially which is due to the 


reduced load factors of the conventional system. 


$ 





In case c) a fundamental change is observed with respect to case b). 

Base load capacity is significantly reduced and, moreover, coal-fired 
plants now replace the greater part of nuclear capacity. This change 

leads to reduced fixed costs expenditures and increased fuel costs, 


whilst the total costs per produced kWh drop as compared with case b). 


Table 9 shows that a simple minimization of total costs may result in 
undesirable fuel substitutions. For the optimized system of case c) 

the larger part of the fuel savings (with respect to case a)) are ob- 
tained for coal and uranium, whereas the savings for the strategic 
fuels oil and gas are negligibly small. This may be contrary to the 
prevailing national energy policy to reduce the oil/gas share for 
electricity production. Energy policy considerations are es: ential in 
realistic system optimization studies. 

For the re-optimized plant mix (case c) in table 9) the concept of 
"capacity credit" looses its significance. Apart from the true capacity 
savings there is a substitution of large units, which contribute re- 
latively little to total system load carrying capability, by smaller 
ones. The apparent capacity credit of wind power is much larger than 
the value found in chapter 3, where a simple one step replacement was 
assumed. In this case the value cof wind power can no longer be assessed 
by separately considering capacity savings and fuel savings as was 
implicitly assumed in chapter 3. 

If the examples of the optimized cases a) and c) are considered as two 
options for a far away future, i.e. a system with and a system without 
wind power, the value of wind energy can be derived by comparison of 


the cost figures in the two cases. Total costs savings are 1029.4 - 


904.8 = 124.6 x 102 Dfl/hr; wind energy production is 1.454 GWh/hr; 
454 20.086 Df1/kWh. 


the break-even costs of vind energy are 124.6 x 1073 / 1. 
Table 9 shows that in this particular ca:e the savings in capital ex- 
penditures substantially contribute to the break-even costs of wind 
energy. 

Break-even costs were calculated in this way for various levels of 
wind power penetration. Results are given in fig. l4, together with 
total costs per kWh produced by the conventional units (e.g. 14.75 ct/ 


kWh for the optimized system with 4000 MW wind power, see table 9). 


LI 


Calculational results for optimized power systems consisting of plants as specified in table 8. 


Table 9. 





Mix optimized with 4000 MW vi 


Varia 


103/hr 


Dii. 


904 .8 





£1. 
03/hr 


376.3 





case c, with 4C 





duction kosts 
GWh/hr ի 
| 


6.134 





103/hr 


Dfl. 


528.5 





Installed | Fixed | Pro- 
capacity| costs 


TT. a 





Plant mix optimized without wind power 


case b, with 4000 MW w 


nd 


* 
` 
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costs 
Dfl. 


103/hr 


465.4 


116.9 


928 .8 





Variable| Total 


103/hr 


165.4 


43.4 


307 .2 





duction lcosts 


Pro- 


GWh/hr IDf 1 . 
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without wind 
Variable| Total 


costs 
Dfl. 


103/hr 


479.6 | 4. 


145.6| 0.543 


1029.4] 6.134 
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72.1 
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The relative contribution of capital savings to the break-even costs 
appeared to increase with increasing penetration level, which is 
almost contrary to what one would expect if capacity credit and fuel 


savings are considered separately. 


4.3. Conclusions 





The LOLP method seems to be quite useful to investigate the load of 
the individual plants and total fuel consumption. The consumed fuel 
mix and possible savings due to wind power integration are dependent 
on the adopted lcad ordering strategy. This strategy may depend on 
the fluctuations of the wind potential because of control limitations 
in the operation of the conventional plants. This suggests that cal- 
culations with simulation models, mentioned in chapter 1, in which 
the control of the conventional plants is modelled in detail, are 
necessary to verify the results obtained with the LOLP method as to 
the feasibility of the chosen strategy. On the other hand, the LOLP 
method may provide useful information on input parameters for simula- 
tion models. Generally speaking, probabilistic models may be less 
accurate than simulation models, but they are more versatile and 


cheaper (with respect to computer costs). 


Results of costs optimization studies show the same trend as found by 
the simulation models mentioned earlier. From the point of view of 


costs optimization as well as stability control of the grid, the in- 


tegration of wind power beyond a certain capacity requires the replace- 


ment of large base load plants bv smaller and better controllable 
plants. Owing to this replacement the concept of displaced capacity, 
discussed in chapter 3, can no longer 5e applied for assessment of 
the value of wind energy. 

Costs optimization leads to undesirable fuel substitutions and should 
be pursued with prudence. The development of small versatile coal 
fired plants could facilitate policies towards fuel diversification 
when a sizable amount of wind power would be integrated into the na- 


tional power supply system. 
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5. SUMMARY 





Studies concerning the integration of wind power into the utility grid 


can be performed by simulation models as well as by probabilistic 
methods. In this report probabilistic methods have been applied to 
estimate the capacity credit of wind power for a single replacement 


mode. 


Constraints with regard to cycling of the conventional plants were 
neglected. An ideal control system was assumed to exist such that 
technically available plants are set on line immediately whenever 
they are required. This assumption may lead to an overestimation of 
the capacity credit of wind power. The degree of overestimation will 
depend on the control properties of the constituent conventional 


plants and the applied control strategy. 


Within these limitations chapter 3 deals with several methods to 
calculate the capacity credit of wind power integrated with a given 


conventional supply system. 


One method applied in chapter 3, the LOLP method, needs special 
mention. The changes in the production patterns of the individual 
conventional plants, due to the contribution of wind power, can be 
studied by this method, subject to the adopted strategy by which 
wind energy is to be absorbed by the grid. It is even possible with 
this method to optimize the power plant mix for varying levels of 
wind power penetration. In that case the concept of capacity credit 


looses its significance. 


It is our opinion that the integration of a sizable amount of wind 
power involves more than fuel savings and displaced capacity alone. 


The changing production patterns of the individual plants and pos- 


sibly the changing structure of the power system are important factors 


for the assessment of wind energy as well. 
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APPENDIX.  Loss-of-load met 





Ճ.1. A load duration funct: 
notes the time fracti: 
demand (or production; 
obtained from the prol 
load demand (or produ: 

F(2) = | p(x)dx. 
Starting with the den: 


duration function of | 


It 1S common practice 


load 
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prob. density p(:) 


1 1 
L Oad 


` 














time fraction F(?) 


-hod. The area below the loac duration curve represents the total 
energv required during the period under investigation. If 
ion, F(2), of demand (or production) de- plants were available all the time, a plant with capacity 
yn of a period (usually one year) with | could produce the energy represented by the area below the 
| greater or equal 2. This function is | dotted line of fig. 1(b) . The load duration curve of this 
ability density function (2) of the plant is shown in fig. l(c); the du-ation function of the 
'tion): TE remaining demand is shown in fig. i(d). A secoud plant can 


now be scheduled and the above procedure repeated. 


This method can no longer be applied if the full capacity c of 
sity function of fig. l(a), the load i the plant is not available continuously. In tnat case one can 
19. l(b) is obtained by integration. define a probability density py (X) for the capacity on outage. 
to plot the load along the vertical axis. The method of section A.l. can now be generalized by considering 

the duration function of the sum of the load % and the capacity 


on outage o. This sum is called "effective load" e: 


f 2 and o are statistically independent the d ration function 
of the effective load, F(x), is the convolution product of 


po (X) and F(x): 


remaining load | (d) 


ーー ニー ニー Շ 


ーー 


F 。(x) f p. (x')F(x-x')dx'. 
օ 


This operation is shown in figs. 2(a) through 2(c). 











m It can be proven that the function F。(x) above the dotted line 
| T - in fig. 2(d) is the duration function of the remaining load 
p" cd U l | after the plant with capacity c and outage probability pg(x) 
-— "hu supplied power i has been scheduled. Fig. 2(e) is thus the generalization of 
E fig. l(d). The load duration function of the plant (cf. fig. 
1 (c) | l(c)) now becomes the product of the functions F(x) and F,(x), 
Ip where F,(x) represents the duration function of available ca- 
pacity (c-0): 
—P Դ 








time fraction eer C C-X 
Fa(x) = f p, (c-x')dx' = | p(x" dx’. 
x Չ 
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T" | however, are either available with probability p or on outage 
The time fraction with remaining demand > 0 (see fig. 2(e)) 
| i24 «Ա with probability զ. In this case we may write: 
is called the loss-of-load probability (LOLP); it is equal 


to F。(c) . The area under the remaining curve, when multiplied Po(x) = p 5(x) *q5(x-c) with p+q Հ 1. 


fte) by the duration of the period under investigation, represents 


- " sun բո " 
atkhobive inka à the energy not supplied ("loss-of-energy"). ատ նձ: 
All generating units can be added one-by-one by successively | Բ (x) p F(x) +q F(x-c) 
e «* し プ . 
repeating the procedure described in section A.2. Each time 
capacity on 


| outage x | a new unit is added by convolving the effective load duration [n this case the load duration function of the plant is equal 
| à 


function with the outage probability density of the new unit, to p.F(x) for Osxsc (and zero outside this range). シン 


which results in a new effective load duration function. For the integration of wind power, however, the general formula 
Finally, a LOLP and a "loss-of-energy" remain, which are the hás to be used. 7 0 


reliability indices of the supply system. 
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The convolution operations have commutative as well as associa- ~t 





Ս time fraction 1502) density ք (x) 0 time fraction IF (e) 
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tive properties and the reliability index LOLP is therefore vt / 
independent of the order in which plants go on-load. This is | 6 
not so, however, for the load duration functions of the indi- 

vidual plants, these functions change by altering the loading \ 


Order 。 


In the strategy in which plants are ordered, it ay be desirable Հ "- 


to partially load a plant initially at, say, 50% of its maximum es | Հ ua, aai, 


』 remaining load | capacity, and to save the remaining capacity for a higher level i bt r 
"^ 


max Vn in the loading order. The method can easily oe applied to adopt 
such a strategy. When a plant is convolved into the system for 
the first time and for only a part of its nominal capacity, the 


outage probability density function for this part of its nominal 


on capacity is used. When the remaining capacity of the plant is 


0 LOLP — 


time fraction system by deconvolution and then re-convolved with its full 











required at a higher level, the plant is first removed from the 











| ' ; i ; i ' 。 d ! 
capacity (this manipulation is allowed by the commutat1Ve and | un -f 
© © 


0.3- 
0.2 - 


time fraction F, 
ե. 
associative properties of the convolution product). This idea 


has been applied in section 4.1 (variant 2). 


If the outage probability density function pg(x) of a plant 
covers the whole range of capacity values from 0 to nominal 
capacity c, the convolution product given in section A.2. has 


to be calculated by numerical integration. Thermal plants, 
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Fig. 13. Production shares of the different plant types of table 6 
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Abstract |. INTRODUCTION | Various operating restrictions on thermal plants cannot be taken into 


account, therefore. 
Է x4 bui | | Integration of wind power into the electricity grid affects total 
ն ower into the electricity grid affects total i | | 
nm M— din The frequency and duration (FD) method can 5e considered as an ex- 


consumption of the thermal supply system as well as total installed fuel consumption of the thermal supply system and the economically 


tension of the aforementioned LOLP method in that it retains some 


; ; : imal mix of primary fuels consumed. It also affects total in- 
る ՉԻ - r a reliable power supply. Both ef- optima ) 
A Frequency and Duration Method for the production capacity necessary fo P pply details of the time-dependent behaviour of system load and production 





ad | | tei : ilisti r model based stalled production capacity necessary for a reliable power supply 
Evaluation of Wind Integration. fects can be quantified with a probabilistic compute 





units. The theory of the FD method for reliability studies of all- 
d the optimal mix of plant types. The effects on fuel consumption 
duration method. -—-— : hà há | 
on the frequency and thermal systems can be found in literature in various formulations; 
and generating capacity are mutually dependent. 


A.J. Janssen The theory of the method, applied to wind-assisted systems, is out- 


see, e.g., |20-22|. In |13| and |16| the FD method has been applied 
る TM have been published dealing with these 
i 1 litv can be expressed in various reliabilit In recent years many papers p g | ՓՈՐԻ 
Netherlands Energy Research Foundation ECN — — — ա " in different approximations to the reliability evaluation of wind- 


indices, e.g. loss-of-load probability and frequency of loss-of- effects of wind power integration. Hour-by-hour simulation models Mais 00500 haniai յկձեննեն) նեխած Libia WIN EAA 
load events. Calculation results show that the different indices have been used mainly to study the effects on fuel consumption with of the additional information offered by this method with respect 
are affected differently by wind power integration. Thermal capacity dua WeeuM HMM OF Th CGS անառ vases voti to the frequency and duration of loss-of-load events. 
savings due to wind power are therefore somewhat dependent on the Probabilistic models have been used mainly to get information on iis cones iube ak 本 seed dep seliabilitp aad MM) լաս» 


adopted reliability criterion power supply reliability of wind-assisted systems and on the capa- 


duction studies of wind-assisted systems. In section 2 the theory 


à; ; city credit of wind power |3-5,10-19|. Most of these statistical 
: am ween wind ener անը, ' ] l . 
Furthermore, there is no unambiguous relation betwe 57 is described; in section 3 some results of calculations are presented, 


production and the type of fuel that is saved. It is shown that even models consider the (annual) loss-of-load probability (LOLP) as the 


which demonstrate the additional information offered by the FD method. 


, ՛ ՛ ՛ - ՛ - < で で (^ ՛ 2 r 1 , 
for high wind power penetration levels a reasonable freedom exists reliability index of the supply system. The LOLP is evaluated by 


in thermal planı operating strategy to achieve savings in accordance combining the load distribution function and the probability density 
1 し 


with the prevailing fuel mix consumption policy. function of available capacity (i.e. capacity that is not on forced 
or scheduled outage). Wind power is treated either as a negative 
load (subtracted from the system load on an hour-by-hour basis) or 
as a multi-state production unit with many derated states. Some of 
these nodels have been used also for probabilistic production 
costing calculations and/or for evaluation of ti re-optimized mix 
of plant types 15-19... In these applications the statistical models 
offer the advantage over simulation models of being easy to handle 
efficient with respect to computer use. However, in these 
applications chronology is lost, since use is made of the load dis- 


tribution function instead of the chronological load function. 


- 


՛ 
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2. THEORY 
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2.1. Basic principles 





The input to the FD model calculations consists of the frequencies 
and durations of system load levels and forced outages of the pro- 
duction units. For a predefined merit orde-, the units are 
scheduled one by one to serve part of the loac, starting with base 
load and ending with peak load. The prime output of the model 
calculations consists of the frequencies and durations of the load 
levels of the various units and of the load that cannot be served. 
In the process of unit scheduling the so-called effective load is 
the key variable. It is the sum of system load and capacity on 
forced outage of the units already scheduled 1$,23,24|. A unit 

can be conceived to consist of a fictitious perfectly reliable 
unit and a fictitious random load with statistical properties equal 
to those of the capacity on forced outage of the real unit. If 

this unit is scheduled its nett production will be the difference 
between the fictitious production of the perfectly reliable unit 
and che fictitious load demand. This load demand is added to the 
system load (which results in the effective load) and may be served 
by units to be scheduled next. Finally, after all units have been 
scheduled some unserved load will remain which gives information on 


տ. 2 $ e4 * 
system reliaDlililty. 


In what follows the algorithm for the evaluaticn of tne effective 


ーー 
E 
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load will be given. First the simple case will be considered that 
conventional production units are scheduled, which are modelled as 


1 


two-state Markov processes; fig. l depicts the Markov process. 


Next the generalization to include multi-state wind power plants will 


be outlined. 









nological load data files will not be accounted for. 





where 

















i 3. Y data for the time period under in- T も i | THF | LE 
Chronological (e.g. hourly) load dat : In studies for expansion planning se will have to be made of forecast 


Po(x) = probability of capacity on outage equal to x, 


























i i . From this data set two functions are bági Է" ՊԱՊԸ ' 
vestigation must be available. From chis M load data which are inherently un.ertain. Dutch utilities assume, in 


q = unit forced outage rate (see fig. 1), 
















- capacity of the unit. 







derived, viz.: 


(2 


their planning studies, an uncertainty of 3.22 in national annual peak 
























1 1 | Fa ( | obability =á | | 1 | 
|) The load distribution function F(E), denoting the probability load. This uncertainty can be taken into account in the load functions 







The independence of i and © also means that the probability of simul- 













F;(E) and £;(E). If one assumes that all forecast load levels of the 
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EN 


F; of load levels greater or equal E. 








taneous transitions of l and o is negligibly small. Hence, the 
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period under investigation have identical fully correlated uncertainties, 




















2) The transition frequency function f,(E), denoting the average 
transition frequency f,(E) is made up of two independent contri- 














frequency f; that levelE is crossed 18 an upward movement of the load. which are normally distributed with standard deviation 7 equal to a 


butions. the first one due to transitions of 2, viz. 
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The two functions are derived by simple counting techniques. For the certain fraction of system peak load, this uncertainty 15 accounted 









daily recurrent load pattern of fig. 2a it is easily seen that the ollowing simple convolutions: (1-q) £, (E) + q £,(E-O), 








corresponding functions F (E) and f, (E) are those of fig. 2b and fig. 
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and the second one due to transitions of o, viz. 









2c, respectively. The area under the distribution function, when 











E-C 






] e 3 
à F(x) = = {F (E-C) -F。(E)」 >» 


multiplied by the duration AT of the period under investigation, de- 
T Լ P 


t " ե n 

£,(E) = | N(x) f£, (E-x)dx, (3) 1 լ 
T し 

E 







notes total energy demand: 
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Վ (x) being the Gaussian probability density function witn zero mean ' à i à 
p" 5 ՛ " mum -- ան, where T is the unit cycle time (see fig. 1). 


































‘max and standard deviation In practice this uncertainty hnas a noticeabie 
En = AT] F,(E)dE, - (1) Thus: 
0 0 A - | | " "un à; Պո 
ZEND .. զա. Օ ։ 0 0 Օ ՕՁ uu աա նն ն | 7 m ան 2 £ (E) = (1-9)£,(E) + q£,(E-C) + (F, (E-C) - F, (Œ)}/T. (5) 
と 









Fe(E) and E。(E) can be evaluated in a recursive way for decreasing 







where Lmax is the system peak load. 


values of E by using the boundary conditions 





As stated earlier, wind power can te modelled either as a negative 















has been scheduled to serve, as far as possible, 





load or as a multi-state production unit. The latter modelling will 
0 for E greater than the maximum load level Lmax» 


F,(E) = £4(E) 











be discussed later. In the former case the differences between syn- 
F(E) = 1 for E <0, 












chronous values of load and wind power should be used in the com- 










- 上 E<0Q. 










putation of the functions F;(E) and f,(E). 






After the unit has been scheduled the remaining energy demand is 
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is shown in fiz. 1. Since transitions of 2 and o occur independentiy, 
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- In computer applications discrete values of the variable E must 













be used. The step-size ^E is a compromise between desired accuracy E) of the effective load can be derived 
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and calculation efficiency. Load fluctuations smaller than AE are 
















neglected because of this .umerical approximation. Also load fluc- | -— A 
5 PP The (expected) energy delivered by the unit is Eg-E, (see (1)). 






The next unit to be scheduled is needed during a fraction of time 







tuations on a time scale smaller than the time-step in the chro- 


equal to F,(C), the (average) frequency with which this unit is 





called for energy production is f,(C), and the (average) duration 


of the periods dur: 
By substituting Fa 
load functions bel 
can now be evaluat: 
scheduled successi 
indices can be deti 
1) Loss-of-1load pr: 

LOLP = F。(Cr 


2) Average frequen: 


ni 


3) Average duratio! 


4) Expected energy 


Average power d 





there is no 
mittent power sour 
multi-state 


characterl 
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) The wind power 
noting the prob 
of wind) equal 


The frequency 
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է, 


that wind power 


ing which this unit is needed is F,(C)/f£,(C). 
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(E) and Fa。(E) for F,(E) and f, (E) in (4) and (5) the 
nging to the effective load with two units scheduled 
d. After all units of the supply system have been 
zely (total capacity Cp) the following reliability 
ormined: 

obability: 


). (7) 


cy of loss-of-load events: 


(8) 
n of loss-of-!oad events: 

(9) 
deficit 
nax 

F_(E)dE.。 (10) 

eficit during loss-of-load events: 
ՆՔ). (11) 





relation in time between system load and an inter- 
ce like wind, the generalization of (4) and (5) for 
tors is in fact quite straightforward. We define 
functions of available wind power. 

outage probability density function p,'/E), de- 

ability p, of wind power on outage (due to lack 

unction £,(E), denoting the average frequency է 


outage levei E is surpassed. 


et 
(n1 
いり 





rm a TI A a ameters. 





These two fuactions are derived from chronological (e.g. hourly) 
wind power data in essentially the same way as the system load 
functions discussed earlier (with the same step-size 3E). Wind power 
fluctuations smaller than AE are thus neglected, as well as fluctua- 
tions on 3 time scale smaller than the time-step of the chronological 
wind data files. 

Making use of the independency of system load fluctuations and wind 
power fluctuations, we now Obtain the effective load distribution 


function with the convolution (cf. (4)): 
F (E) = | p。(x) F,(E-x)dx, 


and the effective load transition frequency function (a direct 


generalization of (5)): 


fe(E) = f p, GO Ex(E-x)dx + | £4(E-x) dF, (x). (13) 


The calculation procedure used to evaluate (12) and (13) is basic- 
ally the same as the one used for (4) and (5), although the cal- 
culation effort is much »igger. 

The assumption made in (12) and (13) with respect to transition 
independence of system load and wind power will be discussed in 


section 3. 


)perat1ng considerations 
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Loading of production untts 





The above described unit commitment policy is not quite realistic. 
In the daily or weekly production planning account is taken of the 


unit conversion efficiency as a function of (partial) output power, 


—~ 


of unit ԵՏ ԱՀԱ: of spinning reserve requirements, etc. AS a 
consequence it may be that during hours of minimum load several units 
will be run at their technical minimum, thus avoiding start and stop 
procedures for these units. Such operating policies naturally in- 
fluence the energy productions of the various units. Partial loading 
of the units can be modelled by segmentation of the unit output power 
range, for example: 

segment | from 0% to 40» of rated power, used to serve base load 

during hours of minimum demand; 

segment 2 from 407 to 90% of raced power, used for daily load 

following; 

segment 3 from 902 tc 100» of rated power, to be considered as 


spinning ceserve and used only in emergency situations. 


The evaluation of the effective load properties with segmented commit- 


ments is a bit complicated since the random outages of different seg- 
ments of a unit are correlated. The easiest way to take these cor- 
relations into account is as follows |23|. Suppose i segments of a 
unit (total capacity C;) have been committed. The effective load is 
described by Ffi(E) and £*1(g). Segment i*l is to be committed. 
sgments are first deconvolve!, which re- 


functions F。+(E) and £5. (E): 


= "Ll 
- ( 
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The functions . 


E) and fo (E) can be obtained in a recursive way 


for rising values of E, making use of the relations 
F'i(g) = 1 for E<0 


£$i(E) = 0 for E<0 


Fti(z) = ԼԵԼ) = 0 for E>Lnax+Cc 


C+ being total capacity already committed. 

In production costing calculations the conversion efficiency depen- 
dence on (partial) output power can be taken into account by defining 
heat rates appropriate to each segment. Furthermore, since for a 
given commitment order the number of starts of a unit can be derived. 
from the effective load function Fa。(E) , total start-up costs may be 
determined as well. 

Finally, it may be mentioned that althougn the adopted commitment 
rule affects the energy productions of the various units, it has no 


effect on the ultimate system reliability indices. 


Start-up frequency Limitations 





Frequent starts and stops of base-load units may be disadvantageous 
economically or may he even impossible technically. One may there- 
fore set a limit to the frequency of starts of such units. If, for 
a predetermined order of commitment, a unit would have to be 

started too often it should not be committed at all during the time 
period under consideration, but it should be replaced by units with 


more flexible operating characteristics. 


The relevant start frequencies are contained in the load function £,(E). 


Commitment priority of wind power 





Since O&M costs of wind turbines are very low it seems logical te 
place wind at the top of the merit order. However, in the case of 


high wind power production during hours of low system demand this 


policy may lead to the necessity to stop large base-load units. 
Because of the limited operating flexibility of these units the re- 
sulting fuel consumption may not be an economic optimum. Therefore, 
it may be better not to place wind (modelled as a multi-state pro- 
duction unit) at the top of the merit order, but first commit a 
certain amount Lmin of thermal capacity. If wind is considered as 

a negative load and is subtracted from the system load one can set 
a lower limit L4;4 to the nett load. The functions F。(E) and £, (E) 
belonging to this nett load then have to be corrected for E Հ Լր 


as tollows: 


Fig. 3 shows this correction 
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3. CASE STUDIES 


The theory described in section 2 has been implemented 
program. Some results of calculations will now be given t 
the various ways in which themodel can be used. In sect 
relicbility indices of wind-assisted systems will be g: 
section 3.2. fuel substitutions due to wind energy pro 


be considered. 


3.1. Reliability indices 





Reliability indices have been calculated for a hypothe 
system and a load demand that have been used and descr: 
in an earlier publication |18|. We will give here only 
of the system and refer to |18| for further details. 
The parameters of the conventional generating system a 
table 1. In order to obtain information on the relatio 
reliability and cotal installed capacity, the number o 
was varied in the calculations. Scheduled maintenance 
sidered. 

As to system load, the half-hourly load values realize 
Netherlands in 1978 were scaled up by 24% to values ex 
with system peak demand 11472 MW and annual average de 
An uncertainty of 3.2% (1 c) in the forecast peak load 
account in the way described in section 2.1. The step- 
the numerical calculations was 25 MW. 

Hourly wind speed values of Den Helder weather station 


1965 were used, together with the power curve of a sta 


turbine |18|. Some properties of the turbine are given 
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in a computer 
.O demonstrate 
ion 3.1 

iven and in 


luction will 


trical generating 
ibed in detail 


the essentials 


re given in 
դ between system 
£ 50 MW units 


was not con- 


d in the 
pected in 1985, 
mand 7588 MW. 
was taken into 


size used in 


in the year 
ndard wind 


| 1n table 2. 


Conventional generating system parameters. 





Unit 


| 


capacity 


(MW) 








Forced 


outage 


rate q 


Mean time 
to repair 


(days) 





nuclear 
coal 

oil/gas 
oil/gas 


oil/gas 





variable 














Wind turbine characteristics. 





| Hub height 


| Cut-in velocity v 


| Rotor diameter 


| Rated power 
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| Lut-out velocity Veo 
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Average wind speed (1965 data) at hub height was nearly 9 ms~! 


(Den Helder is the most favourable Dutch coastal wind turbine site), 
the annual wind turbine load factor was 0.38, the time fractions 
with zero power and rated power output were 0.24 and 0.19, respec- 
tively. In most of the calculations described below wind was con- 


sidered as a negative load. Any effect of dispersed siting was 


neglected. Technical outages of the turbines were not taken into account. 


Base case 





Fig. 4 shows the five reliability indices, defined in section 2.1, 
eqs. (7)-(11), calculated for the thermal system alone and for the 
thermal system assisted by 5000 turbines. It can be seen that the 
system LOLP, the frequency of 1oss-oE-1oad events f and the energy 
deficit Eg are strongly dependent on total installed capacity. 

The average duration of loss-of-load events T and the average power 
deficit P are much less dependent and are therefore less useful as 
measures of system reliability. Figure 4 also shows that the integ- 
ration of wind power affects the different indices in different ways. 
En is reduced less than LOLP, hence the somewhat higher value of P. 
The reduction of f is significantly less,which is reflected in the 
decrease of T: wind power output fluctuations clearly occur more 
frequently than system load fluctuations or thermal unit forced out- 
ages. Because of these differences the capacity credit that may be 
assigned to wind power depends on the adopted reliability measure 
(the capacity credit is the nesteunent distance between corresponding 
curves in fig. 4 at a given level of a reliability index). 

In table 3 some possibilities are given. If LOLP is adopted as the 
reliability measure the capacity credit would be 463 MW, which would 
result in an increase of 11Z of the frequency of loss-of-load events. 


If this frequency would be the adopted reliability measure, the capacity 


Table 3. Reliability indices of a generating system assisted by 
5000 MW wind power, for various values of the capacity credit 
assigned to the wind power plant. The capacity credits refer to 
50 MW units with SZ forced outage rate. Without wind assistance 


the systea contains 81 units of 50 MW (see fig. 4). 





Without 
5000 MW wind power 
wind 








Capacity 











credit (MW) 
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credit would be 50 MW less. This reduction of 50 MW may still be too 
optimistic since wind power fluctuations on a time scale less than 
an hour have been negiected. On the other hand, by geographic dis- 
persion of wind farms local fluctuations on this time scale would be 
smoothed for the greater part. In reality 5000 large wind turbines 
would require an area with dimensions of at least several tens of 
kilometers. 

From the numerical examples given above one may conclude that the 
effect of large-scale wind power integration is such that the use of 
only one index to quantify system reliability oversimplifies the 
problem of system reliability. 

It should be noted that the capacity si-ings have been expressed here 
in terms of small 50 MW units for convenience only. In reality, large- 
scale wind power integration calls for a re-optimization of the con- 
ventional plant mix, which would result in a marked reduction of 
large base-load units and probably some additional peaking plant 117, 
18!. In that case a definition of wind power capacity credit is no 
longer trivial, although the capital savings due to wind power may 


well be indicated |17,18|. 


ar ーー - © Գ. J 
- X と の ん レン Ww 212 





The cut-out wind speed v.g of the turbines has been assumed to be 

20 ms" !. The time fraction with wind speeds higher than Vig is 

1.82, according to the 1965 Den Helder wind data file. If veg could 
be increased to a sufficiently high value, the 5000 MW wind farm 
could thus produce an additional amount of energy aqual to 90 MWy 
per year. The effect of this increase on system reliability is shown 
in fig. 5 (compare the full lines and the broken lines). The in- 
dices LOLP and Eg are affected equally; the reduction of the fre- 


quency f is a bit more pronounced (hence the higher values of T), 
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which can be ascribed to the relatively short durations of periods- 















as follows. If wind power is treated as a negative load,only fluctua- 


was made with assumed mttr values of all units equal to 10 days. 





Wind power and thermal plant data 
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with wind speeds above 20 ms . The capacity credit associated with tions of the nett load are considered on a time scale of one hour 
Wind energy data are again based on the 1965 Den Helder data discussed 


Because the forced outage rates were not altered, the calculated 









or more. If wind is modelled as an independent power source, fluctua- 






the additional energy output is about 50 MW. 










LOLP, Eg and P remained unchanged. The effect on է is shown in fig. in section 3.1. The load factor of the wind turbines during the 




















Lo tions of the nett load within an hour will also occur. In the latter 
“ower ana -ՕՃԶ demand + i " ' ' à ' է 
- | 7 (dotted lines). Since the reduction or 
modelling simultaneous changes of system load and wind power are F 


winter week under consideration was quite high, viz. 0.456. 






for the wind-assisted 





















system is almost the same as for the thermal system alone, and be- Figure 8 shows the system load distribution function F(E) and the 
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Diurnal and seasonal variations of wind availability in the Nether- 
















neglected. Since these simultaneous changes do happen in reality the 


cause this change of f is most probably an overestimation of the transition frequency function f(E) for the week under consideration. 










lands are such that there is a small positive correlation between 



































resulting frequency of nett load fluctuations (and of loss-of-load 


wind power and system load 18| 。 This correlation was taken into 7. - T : ! E E i 
effect of uncertainties in mttr, one may conclude that accurate «now- In fig. 9 these functions are given for the nett load after subtraction 





events) will be an overestimation. One might argue that the two ways 













account in the above examples where wind turbine output power was 1 - bso] isi - E 6 i 
ն ledge օք mttr is not an absolute prerequisite for wind integration of the output of a 5000 MW wind plant. In the latter case medium and 





hidni է T of wind power modelling supply the error bounds of the "true 

ubtracted from system load on an hour-by-hour basis. If, on the | T 

' f : low load levels occur much more frequently, which indicates that ther- 

hib i P | " l l quency f. Although further study would be useful to arrive at a 

other hand, the wind farm is treated as an independent multi-state ai staat epérétioh prob! may be expected (For the sake of 
mal p ems 。 G I 


fundamentally sound definition of the load frequency characteristics, 





power source this correlation is lost. The effect on system relia- ' տ մնա 
| claricy not ali details of the fine-structure in £(E) are shown in 












the error margin shown in fig. 5 indicates that both ways of wind 
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oility is shown in fig. 5 (compare the full lines and the dotted | 
izs. 8 and 9). 











power modelling are quite adequate in system reliability studies. ' | E au d աա 
In this section we consider the fuel consumption of a hypothetical andi r | ' 
The composition օք the conventional generator system 1S given in 







lines). The indices LOLP and Eg are hardly affected, but apparently 












oroduction system during one average winter week in the year 2000. l d に La V | 
| table 4; it strongly differs from the present-day composition and is 


Sustem Load uncertainty 
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the calculated system reliability is improved a bit if the correla- 





Some input data for the calculations were obtained from KEMA (Dutch 





in accordance with a fuel consumption policy far away from gas/oil 






tions are neglected. This unexpected result is caused by the fol- 





Figures 6 and 7 show the influence on system reliability of the as- 

















Electricity Research Laboratory). The data are based on a 





towards coal in large base-load units. AS regards wind power integ- 


lowing effects. At the high level of wi Ow | թ | à; i i 
g à zh level of wind power penetration con sumed uncertainty in system peak load of 3.22. If this uncertainty 




































。 : E scenario with 22 average annual growth of Dutch GNP. Peak load in | - 
would not exist (broken lines in figs. 6 and 7) the same level ot un " , ration this composition can be considered as a worst case. The units 
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sidered here the contribution of wind power to system reliability 


















the average week is 13681 MW, average load is 10238 MW. Ihe pric | | | MM 
listed in table 4 are those expected by KEMA to be available during 


is hardly affected bv av ze win urbi | te 4 : err ; ië 2 , 
" dios: y average wind turbine output but it is mainly reliability would be obtained if installed capacity would be re- 


















っ f heavy fuel oil is assumed to rise 4% per annum to Dfl. 18.50/GJ l | ը | | 
i the week under consideration. We therefore set the forced outage 


determined 5y the probability of こる է | 
bility of non-zero power output (10,18). duced by about 1.42 of system peak load. Neglecting the uncertainty 









rates of the units equal to zero in our calculations. Another reason 


he probability of wind s is above cut-in will incr ՉԱՐ 1 à; | i i 
À speeds above cut-ın wlil increase indeed slightly decreases the wind power capacity credit, for the obvious 
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The present gap between the prices of oil and coal is assumed to 5e 












for neglecting forced outages was to facilitate a comparison of our 
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calculated results with those of a simulation model to be discussed 


wind speeds above cut-out will also increas latt 2 . . ' ; . e 4 i 
: ` also increase. The latter effect integrated with a system which in itself is less uncertain. 













scenario the fu 
later on. 






gative correlation between wind availability and system 



















2 A a na ^ - ーー 1 7 - յ 
て @ ら の の ん ゲ が 211165 72 COMNVENTTOMA L An? TS 


357 of the fuel costs 





The heat rates of the units were approximated by the Willans line: 
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The forced outage rates shown in table 1 are in line with present- 





l f electricity from oil and gas. In order to evaluate the influence 


Although there is hardly any difference between the calculated LOLP 













Q(x) = Q] *Q2.x (16) 






. However, the repair times, 





t ՛ | | 
day experience in the Netherlands |18 4 յ -— n in 
f relative price differences for the dirrerent ԼԱՇԼՏ we also did 


4 m "= 


ւք wind power is treated in one way or the other, 
























mttr, are only reasonable guesses deduced from literature. In order 





prices of gas 






the frequency f or loss-otf-load events does depend on the way in 
| where 















to investigate the effect of uncertainties in mttr, a computer run 





Q = heat rate (GJ/h), 


UJ 


This can be explained 







x = fraction of rated power. 












nuclear 
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nat.gas | 
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Table 5.  Segme! 





| Unit type 
| (rated po 


| of units) 





«ա. 4+ -—— 


| blast fur 


| 2 | nuclear 


| | 

| 3 | coal 

| | 

| 3 | zas 

| 5 | oil 

| 6 | nuclear 

| 7 | coal 

m 

| 3 | gas 

| | 

| 9 | oil 
10 gas 
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-ional system generator parameters. 


The coefficients Օլ and Q2 are included in table 4. They are based, at 



































| | | | ՀՈՐ լ 
Rated | No | զ 1) | 07Ս Efficiency least partly, on performance data of Dutch utilities. The relation 
ja | = ՛ | | 4 
| - | | i Tu o " " 7 a à 
power | OI | (GJ/h) | (GJ/h) at ratos between electricity production and fuel consumption of a unit ls, 
(MW) | units | | power (7) ' 
| | | according to (16): 
առ | | ] 
| | | 
600 | 3 | 450 | 6300 32.0 Qw = 168 F(Ce) Q] + E,Q,/C 
| | | 
600 | 17 | 210 | 5040 |. 45 
| | | լ where 
600 | 2 | 210 | 504 41.1 
| | | 
. | っ | / 
500 | 1 | 175 4200 4].1 Qu = fuel consumption (GJ/week), 
| | 
Լ | | . e / | i TP 
425 | Է | 212.5 3995 36.4 F(C,) = fraction of time during which the unit is needed (see 
| | 
325 | | | 130 3006 37.3 section 2.1), 
100 | 19 | 55 955 35.6 | 
| | կ | լ E = energy delivered (MWh/week), 
100 | 14 | 55 955 | 35.6 
50 | 12 | 107.5 517.5 | 28.8 C - = unit rated capacity (MW). 
| | 
- | | | Start-up fuel costs are neglected in (17). 
explanation of Qj and Գշ. The segmented commitment priority list adopted in our calculations 
is shown in table 5. Blast furnace/coke gas (from the Dutch steel 
nted commitment priority list. works) has to be accepted at all times. The units with rated power 





— | 


greater or equal 500 MW were considered as base-load units and had 




















| Total committed | 
| to be run almost continuously, at a power level greater than or equal 
wer/maximum number | capacity per unit 
| (MW) to their technical minimum (assumed to be 602 of rated power for the 
| ^ 
— 1 nuclear units and 502 for the other ones). These units were not al- 
n | ne 
nace/coke gas (425/1) | 425 | 
| lowed to shut down more frequently than one to two times per week, 
(600/3) | 350 
(600/17) | 300 | else they were not committed at all. No bounds were set to the 
| , 
600/2) | 300 number of starts of the smaller units. The units were allowed to 
(500/1) | 250 "- i ; 
| | produce full power when necessary; the spinning reserve margin during 
(600/3) | 600 
600/17) | 600 | ac peak load hours was very small, therefore. 
600/2) | 600 
(50071 | 500 a Results and dtcusston 
| - 
(325/1) | 325 The rated power of the wind plant integrated with the conventional 
100/19) | 100 | mU | 
| աք | system was changed stepwise in the calculations from zero to 5000 MW. 
100/14) | i | 
| | 
- | - | 
ines (50/12) | 50 | 
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Table 6. Fuel consumption and electricity production during an average 
winter week. 
Electricity demand: 1720 GWh, total fuel consumption without 


wind integration: 16025 17. 








Rated wind power (MW) 





ーー 一 
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| 
| 
| 
| 
| 
| 
| 
| 
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Used wind energy (GWh) 
| Number of committed coal plants 


| Contribution to electr.prod.(2): 


| Fuel costs conv. 


Minimum nett load Lmin (MW) 








Available wind energy EWay (GWh) 


wind 

nuclear 

coal 

zas 

oil 

zas turbines 
consumption (TJ): 
nuclear 


coal 


turbines 


savings (CJ per MWh of 
EW ay) 


use conv. production 





^v | 
(GJ / MWh) 








cost savings ヽ 
『 

1.per MWh of EWay) | 
fee. 


prod. | | 





132.9 

















(D£1/MWh) | 148. 149.7 
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The minimum nett load level after subtraction of wind power (Lmin; 
see section 2.2) was varied also. In table 6 the most important results 
are presented; in figs. 10-13 some results are shown graphically. 

The results (cost savings, fuel substitutions, etc.) are strongly 
dependent on the number of base-load units that took part in the 
energy production. Without wind integration all nuclear and coal 
units were used, as well as one 600 MW gas unit. With 1000 MW vind 
power assistance this gas unit was skipped in the commitment procedure, 
hence the substantial gas savings for this level cf wind power 
penetration. With higher penetration levels several coal units had 
to be shut down as well (see table 6), hence the large coal savings 
for 2500 MW wind. To compensate for this reduction of coal use more 
shall waite had to take part in the production which resulted in a 
nett increase of gas and oil consumption compared with the 1000 MW 
wind power case. À further increase of installed wind power goes with 
a further shift from coal to gas and oil consumption, at least if 
fuel costs (scenario | cost data) are to be minimized. However, from 
figs. 12 and 13 one can deduce that a substantial shift from gas/oil 
consumption to coal consumption can be established at the expense of 
a relatively small amount of used wind energy. At 2500 MW wind power 
penetration this is already of importance. 

The operating penalty of wind power integration is determined to 
a large extent by the shift from efficient base-load units to peak 
load units. The total penalty (including the more frequent partial 
loading of the units and the unused wind energy) can be derived from 
table 6 by comparing the fuel savings (per MWh of available wind 
energy) with the fuel use (per MWh of っ roduced electricity) without 
wind integration. In terms of energy it increases from 6.52 (for 
1000 MW wind power) to about 182 (for 5000 MW wind) of available wind 


energy. In terms of fuel costs (scenario 1 data), however, it is 


23.35 for 1000 MW wind. This negative penalty, or gain, is due to the 
savings of expensive gas. At 5000 MW penetration the cost penalty is 
about 192. The variation of the cost penalty with wind penetration 
level is of course much stronger for the scenario 2 cost data, be- 
cause of the larger difference between gas/oil and coal prices. For 
this scenario the penalty varies from -182 (1000 MW wind) to 237 
(5000 MW wind). 

With respect to cost savings another aspect should 5e considered 
as well. At high penetration levels the thermal plant mix should be 
re-optimized. It is clear that the coal units that could not be 
used for energy production (3 to 6 units of 600 MW for 5000 MW wind 
power, see table 6) should not be built at all, but should be re- 
placed, at least partly, by (cheaper) peak load units (with more 
than 5 coal units shut down the remaining plant was not able to 
serve all the load; in our calculations we assumed that in those 
cases some additional gas turbines were available to serve peak 
load). This shift in investments can lead to substantial capital 
cost savings which are really much more than the capital cost 
savings associated with the capacity credit of wind power discuss«d 


in section 3.1 (17,18 
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Finally, we should like to make a comparison of our calculated 
results with those of a simulation model of the Dutch utilities 
which determines the economically optimal hourly production shares 
of the various thermal units. The model has been run by KEMA with 


the same scenario | input data as we used |26|, however, with more 


detailed specifications of thermal unit regulation range, conversion 


efficiency, start-up costs and minimal shut-down times. The model 








164 


Table 7. Comparison between FD model and simulati 


tions for an average winter week. 





present FD model | sim 





Rated wind power (MW) 0 2500 0 





ーー 


Production snares (ん ) 


w1nd ー 10.6 - 
nuclear ,6.6 16.1 17 
coal 75.6 66.1 76 
gas 7.8 7.4 7 
oil ー 0.0 Հ| 
gas turbines - - «] 
Conv. system efficiency (る 】 38.6 38.3 38 
Total number of start-up 
procedures 17 102 27 
Conv.production costs 
(D£1/MWh) 145.2 145.9 146 
Fuel savings (TJ) - 1578 - 
Cost sevings (Mi 081) - 25.35 - 

















takes wind power into account as a negative load a 
minimum spinning reserve margin of 640 MW. 

In table 7 some figures are given. In view of t 
curacy of our input data and the simple unit commi 
patching rules, the results of the FD model calcul 
prisingly good. With increasing wind power penetr: 
models show the same tendency, but the simulation 
are somewhat more optimistic with respect to, e.g. 
savings. The reason for this discrepancy appeared 
ence in the regulation ranges of the base-load un: 
model a lower limit օք 502 of rated power was assi 


simulation model a technical minimum of 302 was սչ 


on model calcula- 
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ulation model 
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nd maintains a 


he limited ac- 
tment and dis- 
ations are sur- 
ition the two 
model results 
ə fuel cost 

to be a differ- 
ts. In the FD 
med; in the 


sed. This is օք 


course quite attractive with respect to wind integration, although 


utilities would be reluctant to operate their units at such low powers. 


As a consequence, with 2500 MW wind power penetration, only one coal 
unit instead of two was shut down according to the simulation model, 
while all wind energy could still be accepted. This better use of 
base-load units is also reflected in the less dramatic increase with 
increasing wind penetration of start-up procedures of the smaller 
thermal units, although the (too) simple commitment and dispatching 
rules adopted in the FD model will also have contributed to the 
difference in the number of start-up procedures shown in table 7. 


The simulation model results also showed that the number of 


committed base-load units is the key parameter with respect to fuel 


savings and operation penalties due to wind integration, whereas 


the fuel losses associated with the more frequent partial loading 
and starts of the thermal units are of only second-order importance. 
It should be kept in mind, however, that possible penalties 
associated with forecast uncertainty of wind power were not con- 


sidered in the model calculations. 
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4. CONCLUSIONS 


The frequency and duration (FD) method can 5e used for the evaluation 
of power supply reliability and for parametric studies of fuel 
consumption of wind-assisted utility systems. 

Compared with the more frequently used loss-of-load method, the FD 
method supplies useful additional information with regard to system 
reliability. From this information it appears that some caution 
should be exercised in the appreciation of the capacity credit that 
may be assigned to wind power. 


ԵՉ 
However, the two-state model for the thermal production units 


applied in the FD method is still rather simple, the main approximation 


being the assumption that no operating constraints would exist. The 


effects on system reliability of the limited load following capabilities 


and the finite start-up times of the units are thus neglected. 
Recently, a few models have been developed which take into account 
several of these operating constraints |27,28|. Meanwhile, the im- 
provements developed in |27| have been applied by us in reliability 
studies of wind-assisted systems to be reported elsewhere. 

In production costing studies the FD method clearly supplies more 


information than the loss-of-load method concerning the time- 


dependent behaviour of the production units. Some possible restrictions 


to this behaviour can be taken into account in the model calculations. 
Nevertheless, it remains desirable to check the model results by 
comparing them with the results of a simulation model which takes 
all thermal plant operating constraints into account. Such a com- 
parison was made in section 3.2. It was found that the results of the 
two models compare very well, discrepancies could be traced back to 


be caused mainly by differences ir input data. 


+ 


The main advantage of a probabilistic model over a simulation model 

is its very efficient computer use. Therefore, it is very well suited 
for parametric studies. 

The calculated results given in section 3.2. give rise to the important 
conclusion that even for a high level of wind power penetration a 
reasonable freedom exists in thermal plant operating strategy to 
achieve fuel consumption shifts in accordance with the prevailing 
national fuel saving policy. This is in agreement with conclusions of 


Diesendorf and Martin drawn from their simulation model calculations 


ACKNOWLEDGEMENTS 


The wind data used in this paper were obtained from the Royal 
Netherlands Meteorological Institute KMI. System load data were 
supplied by KEMA, as well as the scenario | data used in section 

The discussions with G.H. Bontius, W.A.K. Maas and C.F. Muiser 


of KEMA are gratefully acknowledged. 


REFERENCES 


3. Sørensen, The regulation of an electricity supply system 
including wind energy generators. 2nd Int. Symp. Wind Energy 
Systems, Amsterdan, paper Gl, BHRA Fluid Engineering, Cranfield, 
Bedford, England (October 1978). 

M. Diesendorf and B. Martin, Large Scale Wind Power for Western 
Australia. Proc. Conf. Solar Realities in Western Australia in 
the 1980's, Int. Solar Energy Soc. (ANZ Section) and Univ. of 
Western Australia, p.47-52 (November 1979). 

W.D. March, Requirements Assessment of Wind Power Plants in 
Electric Utility Systems. Final report EPRI-ER-1978, Vol.2, 
General Electric Co., Schenectady, New York (1979). 

Ch.A. Lindley and W.C. Melton, Electric Utility Application of 
Wind Energy Conversion Systems on the Island of Oahu. 

ATR-78 (7593)-2, The Aerospace Corp., El Segundo, California (1979). 
L. Jarass, L. Hoffmann, A. Jarass and G. Obermair, Wind Energy. 
Springer Verlag, Berlin (1981). 

G.E. Whittle, E.A. Bossanyi, C. Maclean, P.D. Dunn, N.H. Lipman 
and P.J. Musgrove, Ճ simulation model of an electricity generating 
system incorporating wind turbine plant. 3rd Int. Symp. Wind 
Energy Systems, Copenhagen, paper L2, BHRA Fluid Engineering, 
Cranfield, Bedford, England (26-29 August, 1980). 

E.A. Bossanyi, G.E. Whittle, N.H. Lipman and P.J. Musgrove, 

Wind turbine response and system integration. 3rd Int. Conf. 


Future Energy Concepts, London, p.296-302, IEE (27-30 Jan. 1981). 


^ X. 





A | E 1 | : 
1 ` Ա 


M. Diesendorf and 3. Martin, Integration of Wind Power into 


Australian Electricity Grids Without Storage: A Computer 


՛ 


t? 


Lon. Wind Engineering, Vol.4, No.4, 1-226 (1980). 

. Whittle, Effects of wind power and pumped storage in an 
electricity generating system. 3rd BWEA Wind Energy Conference, 
Cranfield, Bedford, England, BWEA and IEE (8-10 April 1981). 

Haslett and M. Diesendorf, The capacity credit of wind power: 

չ p.391-401 (1981). 
3. Martin and M. Diesendorf, The capacity credit of wind power: 
a numerical model. 3rd Int. Symp. Wind Energy Systems, Copenhagen, 


paper L3, id Engineering , Cranfield, Bedford, England 


R.G. Deshmukh and R. Ramakumar, Modeling and simulation of WECS 
utility sys 10th annual Pittsburgn Conf. on 
Modeling and Simulation, p.9 
April 1979). 
Deshmukh and R. Ramakumar, Reliability of wind-assisted 
systems, Record of the Control of Power Systems Conference 
79 CH 1377-1 REG.5, p.143-148, 
Texas (March 1979). 
A probabilistic simulation model for the 
calculation of the value of wind energy to electric utilities. 
lst BWEA Wind Energy Workshop, April 1979, Multi-Science 
Co. Ltd., London (1979). 
Fegan and C.D. Percival, Integration of intermittent sources 
Baleriaux-Bocth production Cost Models. Paper A80 111-5, 


“ES Winter Meeting, New York (February 3-8, 1980). 


J.C. VanKuiken, W.A. Buehring, C.C. Huber and K.A. Hub, 


՛ 


Reliability, Energy, and Cost Effects of Wind-Powered Generation 
Integrated with a Conventional Generating System. ANL/AA-17, 
Argonne National Laboratory, Argonne, Illinois (1980). 

M. Diesendorf, B. Martin and J. Carlin, The economic value of 

wind power in electricity grids. Proc. Int. Coll. on Wind Energy, 
Brighton, UK, British Wind Energy Association (27-28 August, 1981). 
A.J. Janssen, T.D. Oei and J.B. Dragt, Statistical methods for 

the assessment of wind power integration into the electricity 
supply system. ECN-101, ECN, Petten, The Netherlands (1981). 

K.F. Schenk, P. Uko and N.S. Rau, Integration of intermittent 
energy sources into a utility resource plan by the method of 
moments. IEEE Communications and Power Conf. Montreal, Quebec, 
p.223-227 (15-17 Oct. 1980). 

A.X. Ayoub and A.D. Patton, À frequency and duration method for 
generating system reliability evaluation. IEEE Trans. on Power 
Apparatus and Systems, Vol. PAS-95, no.6, p.1929-1933 (Nov./Dec. 
1976). 

A.D. Patton, A.K. Ayoub and C. Singh, Power system reliability 
evaluation. Electrical Power and Energy Systems, Vol.l, No.3, 
p.139-150 (October 1979). 

R. Billinton, R.J. Ringlee and A.J. Wood, Power-System Reliability 
Calculations. The MIT Press, Cambridge, Mass. (1978). 

R.L. Sullivan, Power System Planning, McGraw-Hill, New York (1977). 
T.S. Dillon, R.W. Martin and D. Sjelvgren, Stochastic optimization 
and modeling of large hydrothermal systems for long-term regulation. 


Electrical Power and Energy Systems, Vol.2, No.l, p.2-20 (Jan.1980). 


G.R. Fegan and C.D. Percival, Problems in the integration of 
intermittent sources into utility production costing models. 


Int. Symp. on Circuits and Systems, Houston, Texas, 


KEMA Report TEA CM/SMW nr. 80 - 1366 (in Dutch). 

A.D. Patton, C. Singh and M. Sahinoglu, Operating considerations 
in generation reliability modeling - an analytical approach. 
IEEE PES Winter Meeting, New York, Paper A80 082-8 (Febr. 1980). 
K.W. Edwin, H.-D. Kochs and G. Traeder, Untersuchung der Kraft- 
werkreserve in Verbundsystem. Forschungsberichte des Landes 
Nordrhein-Westfalen Nr. 2816/Fachgruppe Bergbau/Energie. 


Westdeutscher Verlag, Opladen (1979). 





F1g.2. 
FIGURE CAPTIONS 
1.1. Two-state Markov model of conventional production units. 

mean time to 

mean time to repair 

mean cycle time = mttf + mttr 

forced outage rate = mttr/T 

1/T failure frequency 

Fig.6. 


System load functions 
a: daily recurrent load pattern; 
b: distribution function g,(E); 


c: transition frequency function f, (CE) 


Correction of the load distribution function F, (E) for the 


lower limit Lin of nett load after subtraction of available 


wind power. ' 
a+b +c = energy demand 


a = accepted wind energy 


Fig.8. 
spilt wind energy 
Reliability indices of a hypothetical power system with and Պր 
Fig.9 
without 5000 MW wind power assistance. LOLP = loss-of-load 
probability; f = average frequency of loss-of-load events; 
Ea = expected energy deficit; T = average duration of loss-of- 
- Fig.lO. 


load events; P = average power deficit during loss-of-load 


events. 


Fig.ll. 
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ibility indices of a power system assisted by 5000 MW 

power. 

—e base case (identical to fig.4) 

ー@ increased wind turbine cut-out velocity 

.@ wind plant considered as an independent multi-state 
power source. 


otes the number of installed 50 MW thermal units. 


bility index LOLP of a power system with and without 
MW wind power. 
—e base case (identical to ք1ջ.4) 


-e system load forecast uncertainty neglected. 


ibility index 1/E of a power system with and without 
MW wind power. 

—e base case (identical to fig. 4) 

-€e system load forecast uncertainty neglected 


.€ increased thermal unit failure and repair times. 


ibution function F(E) and transition frequency function 


of load demand during one average winter week. 


ibution function F(E) and transition frequency function 


rh 
բ) 


of load demand after subtraction of the output ol 


MW wind power plant. 


ction shares of the different production units for 
us levels of wind power penetration, according to 


rio | cost minimization calculataions. 


savings for various levels of wind power penetration, 
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according to scenario 1 cost minimization calculations. 


Savings are expressed in percent of total fuel consumption 


without wind integration. 


Production shares of the supply system including 5000 MW 


wind power foc various levels of minimum nett load Lin’ 


Fuel savings due to 5000 MW wind power integration for 
various levels of minimum nett load. Savings are expressed 
in percent of total fuel use without wind integration. 

: energy savings; -----: fuel cost savings (scenario | 


cost data); .....: fuel cost savings (scenario 2 data). 














avail - 
able 
capacity 





Cr 







































(sADp) aun] -« 


| 




























































system load 








nett load 











5000 MW wind ——— no wir 











ab էի տ. | 
(hr/ y) (year) (MW h/d) (hr) 
| に R 。 
\ 
LOLP \ 
a \ 
T \ 
2.2 ト 
աալ P օ-...ց.... 

















———ÁmÉ- number of 50 





(MW) 





400 


360 








200 





78 80 82 


W thermal units 





















































5000 MW wind di 
で 
. ヒ 
» 
〇 
c 
ի 

58 70 72 7L 76 78 80 82 

— — m» number of SO MW thermal units Y 
で 
Շ 
> 
= 
= 
O 
So 
5 
L の 


























| 









units 












J 
i 
76 
50 MW thermal 


E 
ՎՀ o 
c 
Ս 
O 
Ja 5 
"c 














































f(E) 
1.4 ト 
1.2 ト F(E) 
oh! 
Or 
or 























9 2 L 6 8 10 12 14 
load E (GW) 


























Savings 




















Electr production 9/2 
share % ի Electr. production sa vin 
| à 20 ト share (96) JA 

coal 1 է 

























100 T 
եռ dins _ 。 99$ * oi 


80 


gas + oil 


ーー 






















a 
— sum of 
all fuels 

















—— 











nuclear 











nuclear 










wind 


0 լ լ 
4000 5000 6000 


-m —»- minimum nett load (MW) 


gas + oil 









wind 
2500 5000 
wm rated wind power (MW) 


























0 1000 2500 5000 
= rated wind power 





































0 2 L 6 8 10 
» nett load (GW) 












15 








sum of 

















all fuels 
10 ト E 
iu ^ 
T 
| սազ 
-9F - 
gas + oll 
-10 ト d 
1 լ Վ 
4000 5000 6000 


————m- minimum nett 


load (MW) 








Operating Considerations in Reliability 





Modeling of Wind-Assisted Utility Systems 











A.J. Janssen 











